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Abstract. Discrete spectral transformations of skew orthogonal polynomials are presented. 
From these spectral transformations, it is shown that the corresponding discrete integrable 
systems are derived both in 1 + 1 dimension and in 2 + 1 dimension. Especially in the (2 + 1)- 
dimensional case, the corresponding system can be extended to 2 x 2 matrix form. The 
factorization theorem of the Christoffel kernel for skew orthogonal polynomials in random 
matrix theory is presented as a by-product of these transformations. 
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1 Introduction 

Orthogonal polynomials and their varieties have been investigated for a long time and yielded 
a lot of contributions to various areas including quantum physics, random matrices, wavelets 
and so on. Over the last twenty years, several integrable systems, both in continuous variable 
and in discrete variable, have been found to relate to such polynomials through their spectral 
transformations. One of the well-known results is the relation between ordinary orthogonal 
polynomials (OPs) and Toda chain (see [5, 21] for further details). For recent years, from this 
point of view, new types of integrable systems have been derived [23, 22] and the algebraic 
structures of their solutions also have been discussed in [19, 18]. 

In this paper, we treat the skew orthogonal polynomials (SOPs), first introduced in the theo- 
ry of random matrices [7]. It is known that the continuous spectral transformation of SOPs 
induces the semi-discrete integrable system, known as Pfaff lattice [2, 3], and its relations to the 
Lie algebra, geometry and matrix eigenvalue problems are discussed in [2, 14, 15] respectively. 
Although, the full-discrete integrable system associated with SOPs, namely the discrete coun- 
terpart of Pfaff lattice, still remains to be found and even the discrete spectral transformation of 
SOPs is not obtained. Recalling the application of the discrete integrable system to several areas 
such as eigenvalue problems and so on (e.g. [17, 26]), the discrete integrable system associated 
with SOPs, if obtained, can be expected to have an application or contribution to such areas. 
We shall derive the corresponding full-discrete integrable systems from their discrete spectral 
transformations. 

This paper is organized as follows. In Section 2, a brief review of the relationship between OPs 
and the discrete time Toda equation is given and we propose the discrete spectral transformation 
of SOPs. From these transformations, we derive the corresponding discrete integrable system, 
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which takes (l + l)-dimensional form. In Section 3, another discrete integrable system associated 
with SOPs is shown to be derived also from their discrete spectral transformations, which takes 
(2 + l)-dimensional form. In Section 4, we also show that the system obtained in Section 3 can 
be extended to 2 x 2 matrix form. In Section 5, concluding remarks are presented. 



2 Skew orthogonal polynomials 



We first briefly review a relation between the discrete Toda time equation and the discrete 
spectral transformation of OPs according to [21]. The discrete time Toda equation (on semi- 
infinite lattice) 



^ n* n 
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is known to admit the following discrete Lax pair: 
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for all t. If we normalize 4>q(x) = 1, Favard's theorem shows that {(finix)}^^ is a sequence of 
monic orthogonal polynomials, i.e. there exists some linear functional C 
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such that 
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(x)^(x)] = his* 



Therefore, one can easily find that the discrete Lax pair (1) is just the discrete spectral trans- 
formation of orthogonal polynomials known as Christoffel transformation and Geronimus trans- 
formation [6, 8], respectively. 

Our main aim is to derive the discrete integrable systems associated with SOPs. To that 
end, we shall give the discrete spectral transformations of SOPs like Christoffel and Geronimus 
transformation. We first begin with the definition of SOPs: 



Definition 1. Let (-|-) be the bilinear 2-form from x M[z] to 
symmetric relation 

(f(z)\g(z)) = -(g(z)\f(z)). 



which satisfy the skew 



This skew-symmetric bilinear 2-form is often called "skew-inner product". The SOPs {qn(z)}^ = Q 
with respect to the skew-inner product (-|-) are defined as the set of polynomials satisfying the 
following relations: 



{Q2m{z)\q2n+i{z)) = r n 5 mn , r n / 0, 
{Q2m(z)\q2n(z)) = - (g 2 m+l (z) |<?2n+l (%)) 



(2) 



for all nonnegative integers m, n. 
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From now on, unless specified, we assume that SOPs are monic, which ensures the uniqueness 
of SOPs of even degree. However, SOPs of odd degree have the following ambiguity: the skew 
orthogonality relation (2) is invariant under the replacement 

q2n+i(z) -> <?2n+i0) + a n q 2n (z), (3) 

for any a n £ M. In many cases, the coefficient of z 2n in q^n+lO 2 ) is usually chosen as and then 
the SOPs are uniquely determined. 

The discrete spectral transformation of SOPs (that is, the transformation from SOPs to 
another SOPs) can be constructed. 

Theorem 1. Let {qn(z)}^ =Q be SOPs with respect to the skew-inner product (-|-). Let us denote 
some constant parameter by A satisfying 

A G {z G C | q2n(z) / 0, for all n = 0, 1, . . . }. 

Consider the new sequence of polynomials {<2n( z )}n^=o defined by 

1 sh r n q2k+l{z)q2kW ~ <72fc+l(A)<72fc(2) \ 




<?2n(A) 



(4) 



<?2n+2<» - ^^ly q2n(z) ) + a n ql n {z), (5) 

where a n is an arbitrary parameter. Then {qn(z)}^ = Q are again SOPs with respect to the modified 
skew-inner product 

<.|.>* ==<(*- A) •!(*-*)•>. (6) 

Proof. In this proof, because of the equation (3), we take the arbitrary parameter a n = 
for all n in (5) for simplicity. In order to prove the theorem, it is enough to check the skew 
orthogonal relation (2) for {(/* (z)}^ =0 . Using the skew orthogonality relation (2) for {<7n(z)}^l , 
we find 

<92n(«)>92 m +l( z )>* = ((z ~ A)<?2n ( z ) , {z ~ A)<?2 m+ 1 0=0) 



g2fc(A) g 2 fc+i(^) - q2k+iWq2k(z) 



u , : o r k r n 1 q 2n {\) 
r n / q2k(^)q2k+i(z 



' n \ ^ 
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, , 92m+2(A) , . 

q2m+2\z) T^-q2m\Z) 

q2m{A) 

, \ 92m+2(A) , A ^2n+2(A) 

<72m+2(^) 7T^g2m(2) ) = T^r r nK 

q2m\A) / <?2n(A) 



Similarly, direct calculations also show us 

{q2m{z),q*2n{z)Y = (?2n+l ( z ) , Q2m+1 ( z ))* = °- 

This completes the proof. ■ 

Now observing the discrete spectral transformation of even-degree SOPs (4), one can easily 
find that except for the multiplier factor, this is equivalent to the Christoffel-Darboux kernel 
for skew orthogonal polynomials (we will call this "skew-Christoffel kernel") [9, 16] introduced 
in the theory of random matrices which takes the form 

t i \ Q2k+i(x)q2k{y) - q2k+i{y)q2k{x) , , 

I N (x,y) = 2_^ , (7) 

k=o Tk 

where {q n (z)}^ =0 are some SOPs. We hence obtain the following corollary: 
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Corollary 1. The skew-Christoffel kernel (7) can be factorized into the two SOPs: 

l N (x,y) = ( X -y) qMx)q ™ {x \ (8) 

where the SOPs q 2 N(%) and q2 N (x) are connected to each other by the spectral transformation (2) 
with z = x, X = y. 

There exist many works on the skew-Christoffel kernel for the specific weight cases (see 
e.g. [1, 9]). Here we stress that the factorization (8) holds not only for the case mentioned above 
but also for any well-defined skew inner product. Especially, consider the skew inner product 
in orthogonal and symplectic ensemble given by: 



(/(*)l0(z))orthogonal = / / Sgn(zi - Z 2 ) f (z 1 )g(z 2 )w(zi)w(z 2 )dz 1 dz 2 , 

J Jo 2 

(/(^IsO^symplectic = / (f(z)g (z) - f (z)g(z))w(z)dz, (9) 



where f'{z) = fj{z) and D, w(z) are the corresponding integral domain and weight function 
respectively. For such skew inner product, the modified skew inner product (6) can be written 
as follows: 



((z - X)f(z)\(z - X)g(z)) orthogonal = / / sgn(zi - z 2 ) f (z l )g(z 2 )w(zi)w(z 2 )dzidz 2 , 
{(z - X)f(z)\(z - X)g(z)) symp i ectic = / (f(z)g'(z) - f'(z)g(z))w(z)dz, 



with w(z) = (z — X)w(z). This shows that the skew-Christoffel kernel in random matrix theory 
can be factorized into two even-degree SOPs with respect to the different weight function. 

Remark 1. Similar relations hold between the Christoffel transformation of OPs and the 
Christoffel kernels. Hence we call the discrete spectral transformation (4), (5) as the skew- 
Christoffel transformation. 

Iterating the skew-Christoffel transformation, we obtain a sequence of SOPs by the following 
procedure: 



f n {z) := q n (z), 



n-1 



fc=0 fc=0 



Z — A 

where 

At r i ijkW jjt _*k g 2fc+l( A ) r t 4n+2( X ) 

n ' k r\ qUXY ^ r{- q\ n {X) ' °» q\ n {X) " 

Here we denote the corresponding skew-inner product by (-|-) ■ Then, skew-orthogonality rela- 
tions are given as follows: 

<-l>° — <-l>, M* 44 :=<(*- *) •!(*-*)•>'. 

(4(4 32n+l(*)>* = rt Jmn, {<& m (z) , q\ n (z)) 1 = (ql m+1 (z), q^^z))* = 0. (10) 
For simplicity, we assume at = 0. 
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We then proceed to find the inverse transformation of the skew-Christoffel transformation. 
Using the skew-orthogonality relation (2), we can indeed construct the different contiguous 
relations between {qn( z )}r^=o an d { ( 7n +1 ( z )}^=o as follows: 



n-l 



n-l 



k=0 k=0 

n n— 1 

QL+l( z ) = 4n|l(*0 +Z)7n ) ifc?2fe 1 (^) + Yl 4^2fc+l ( Z ) » 
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Combining (10) and (11), we can form the following Lax pair: 
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where the wave- function = (^(z), g^(z), ^(-z) 



(11) 



(12) 



is the sequence of SOPs. From the 



compatibility condition of (12), the discrete Lax equation is derived: 



L t R t = R t+i L t+i_ 



(13) 



We can observe that if we take the suitable limit, the second equation of (10) goes to the relation: 
d 



dt 



(•!•) = <*•!■> + <•!*■>, 



which appears in the Pfaff lattice proposed by Adler et al. [2]. The discrete Lax equation (13) 
thus can be considered to be the discrete counterpart of the Pfaff lattice. 



3 Discrete integrable systems of SOPs in 2 + 1 dimension 



In the previous section, we derive the discrete integrable system (13) from the discrete spectral 
transformation of SOPs. However the system is non-local and therefore complicated. The 
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similar problem occurs in the case of ordinary biorthogonal polynomials and discrete two- 
dimensional Toda (discrete KP) equation. This problem can be solved by introducing "another 
time-evolution" [25]. Furthermore, observing the relation between skew-Christoffel kernel and 
SOPs (8), one can realize that only even-degree SOPs appear. We hence expect that even- 
degree SOPs play an crucial role in the corresponding discrete integrable systems. Due to 
the guiding principles mentioned above, we shall construct the discrete integrable systems in 
(2 + l)-dimensional form associated with only even-degree SOPs. 
Let {■\-) s ' t be a skew-inner product such that 

{ .\. r +i,t = {{z _ M) . \ {z _ ( .| r ,m = {{z _ A) . \ {z _ A) .)M (14) 

where we assume A 7^ \i. We denote the corresponding SOPs by {^'' (z)}^L . From Theorem 1, 

o are connected by the relation (4) and (5) with A — \i. Similarly, 

{gn*(<z)}£L are also relate( i to {<ln t+l {z)}n=o by (4) and (5). 

We shall find the local contiguous relations among we first recall 

the Pfaffian expression of ordinary SOPs [2], which will play a crucial role in the following 
sections. Pfaffians are defined by: 

Pf(io,...,z 2 n-i) := ^T^" II Pf (V(2i)> V(2i+l)), 

o"£S2n 0<i<n— 1 

where the elements of Pfaffians Pf(i, j) are supposed to satisfy the skew-symmetric relation: 

Pi(i,j) = -Pf (15) 

From the definition of Pfaffians, we can find the following properties of Pfaffians which resemble 
those of determinants: 

2n-2 

Pf(i , • • • , I2n-l) = (- 1 ) k ( i k,i2n-l)(io, ■ ■ ■ , «Jfc, • • ■ , hn-l), (16) 

k=0 

Pf(z , ...,a-i a + b-i b ,..., »2n-i) 

= aPf (i , ...,i a ,..., i 2 n-l) + bPi(i , . i 2n -i), (17) 
Pi(...,s,...,s,...) = 0, (18) 

where j means the deletion of j and 

Pf (t, a ■ i a + b ■ i b ) := oPf (i, i a ) + 6Pf (t, %)■ (19) 

For further details of Pfaffians, see e.g. [13, 20]. 

Employing the "skew" moments (z l \z J ) as the elements of Pfaffians, SOPs can be written in 
terms of Pfaffians as well as ordinary OPs in terms of determinants. 

Proposition 1 ([2]). The SOPs with respect to (•]•) can be expressed in terms of Pfaffians: 
Pf(0,l,...,2n,z) 



q2n{Z) 



Pf(0, l,...,2n-l)' 



Pf(0,l,...,2n-l,2n + l,z) 

q2n+l{z) = q2n{z) = pf( , Q ^ 2w - 1) a ng2n(2), 

where a n is an arbitrary parameter and the elements of Pfaffians are defined by 

Pf(i,j) = (z i \z j ), Pi(i,z) = z i 
for all nonnegative integers i, j. 
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Extending this proposition, we have obtained the Pfaffian expressions of the time-evolved 
even degree SOPs. 

Theorem 2. The even degree SOPs {^W^iq^H*)}^*^*)}^, {?£ M+1 (*)}£* 
can &e written in terms of the moments with respect to ('|-) •' 

, tt Pf(0,l,...,2n,z) 
q2n[Z) Pf(0,l,...,2n-1)' 

s+l,t ( s _L_ Pf(0,l,...,2n+l,/i, Z ) 
92n W z-n' Pf(0,l,...,2n,/x) ' 

S)t+1 _ _J_ Pf(0,l,...,2n + l,A,z) 
l2j z-A - Pf(0,l,...,2n,A) ' 

s+i,t+i, \ = 1 Pf(0,l,...,2n + 2,/x,A,z) 

92ri lZj (z-/i)(z-A) ' Pf(0,l,...,2n + l, / u,A) ' 1 j 

where the elements of Pfaffians are defined by 

Pi(i,z)=z i , Pf(i,A) = A\ Pi(i, / j) = n i , 

Pf (z, n) = Pf (n, A) = Pf (A, z) = 0. (21) 
Before proving this theorem, we show the properties of the Pfaffian elements. 

Lemma 1. Employing the notation (19) in (17), for the elements of the Pfaffians defined 
by (21), the following properties hold: 

Pf ((t + 1) - m • i, U + 1) - » ■ 3) = Pf (* + 1, i + 1) - /xPf (t + 1, j) - /xPf (i, i) + /^ 2 Pf 
Pf((i + l)-/i.»,/i) = 0. 

Proof. Following the definition (19) and anti-symmetric property of Pfaffian elements (15), 
one can find 

Pf ((i + 1) - /i • i, (j + 1) - m • i) = Pf ((« + 1) - V- ■ id + 1) - A«Pf ((» + 1) - M • h j) 
= -Pf (j + 1, (i + 1) - n ■ i, ) + /zPf (j, n(i + l)-fi. i) 

= -Pf (j + 1, % + 1) - M pf(j + 1, i) + M pf(j, » + 1) - /i 2 Pf (j, ») 

= Pf(i + 1, j + 1) - /iPf(i + Id) - fjPf(id) + ^{hj)- 
One can also find 

Pf (i + l-fi'i,n) = Pf (i + 1, fi) - /iPf (», //) = - = 0. 
This completes the proof. ■ 
Proof of Theorem 2. We introduce the function r^' 4 and f£'* as follows: 

<•* = Pf(0,...,2n-l), f^ = Pf(0,...,2n,z). (22) 
It is also convenient to introduce the shift operator T s and Tt satisfying 

T s [f(s)] = f(s + 1), T t [f(t)] = f(t + 1) 
and the notation 
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From Proposition 1 and the relation (14), It is easy to see 



&*(*) = % (23) 

and 

T s [Pf(i,j)} = Pf(i + 1, j + 1) - M Pf(i + 1, j) - M Pf(i,j + 1) + M 2 Pf(i, j), 

T t [Pf(i, j)] = Pf(i + 1, j + 1) - APf(i + 1, j) - APf(i, j + 1) + A 2 Pf(i, j). (24) 

With the help of (17) and (18), we can see 

Pf (0, . . . , 2n, At) = Pf (0, . . . , 2n, /i) - /iPf (0, . . . , 2ra - 1, 2n - 1, /x) 
= Pf(0, ...,2n- l,2n^"l, /x), 

where j denotes (j + 1) — /j • j in the sense of (17). Iterating this transformation, we can find 

Pf (0, . . . , 2n, At) = Pf(0, 0,1,..., 2n^l, fi). (25) 
Meanwhile, we can also find 

Pf(0,i O , . . . ,i2 n -l,[l) 

2n-l ^ 

= Pf (io, . . . , z 2n _i) - ^ (-l) fc Pf(4, /i)Pf (0, ian-i, M), (26) 



fc=0 

where we have used the Pfaffian expansion (16). With the help of Lemma 1 and (24), it is 
straightforward to see 

Pf(7J)=T fl [Pf(i,j)], Pf(i,/i)=0. 
Therefore, combining (25) and (26), we arrive at 

r a [Pf(0,...,2ra-l)] =Pf(0,...,2n,/i). 
In a similar manner, by using Pf (i, A) = (A — /i)A l , we further obtain 

Pf (0, . . . , 2n + 1, /x, A) = A At T s [Pf(0, ...,2n, A)] . 

Finally, arranging these results, we have 

t' +1 «* = Pf (0, . . . , 2n, M ), = Pf (0, . . . , 2n, A), 

A- 1 Pf(0,...,2n + M,A). (27) 



„s+l,t+l 



From these results, we can easily see that f^'* is equal to r^ +1 '* with the substitution z = X. 
Hence, in a similar manner, we can also obtain 

f s n +1 * = z^Pf (0, . . . , 2n + 1, At, z), f *' m = z^Pf (0, . . . , 2n + 1, A, z), 

r s n +W = z'h^X-^m . . . , 2n + 2, At, A, z). (28) 

Finally, from (23), (27) and (28), we obtain the relation (20). This completes the proof. ■ 
Using Theorem 2, we can obtain the following contiguous relations of even-degree SOPs. 
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Theorem 3. The series of even degree SOPs {q2n( z )}ri=0 satisfy the following contiguous rela- 
tions: 

z»,xq£ ht+1 (z) - q 2n+2 (z) = z x C^q s 2 i +1 (z) - z,D^q^\ (29) 

with 



s,t s+l,t+l s,t _a+l,t+l 

As,t _ (,, \ \ 'n+ijn-l r>s,t _ i \\ T ™ T n 

A n -W-A) s+ l, t i^+i ) &n -{fl- A J s+u Sjt+1 , 

s+l,t s.t+1 s,t+l s+l,t 

T n+l T n T n+1 T ™ 

where Tn'* is defined in (22). 

Proof. Using the identities of Pfaffians [13]: 

Pf (io, • • • , *2n-l 5 O, b, C, d)Pf (i , • • • , «2n-l) = Pf(io, • • • , «2n-l, a, b)Pf(i , . . . , i 2 n-i, c, d) 

- Pf (if), I2n-l,a, c)Pf («o, • • • , «2n-l, &, d) 

+ Pf(i , • • • , i2n-i, a, d)Pf (t , . • • , i 2n -i,b, c) (31) 

and 

Pf(i , • • • , i2n, a, c)Pf (i , • • • , «2n, d) = Pf(io, ■ • • , «2n, a, ft, d)Pf (i , • • • , «2n, c) 

- Pf (i , • • • , *2n, a, C, d)Pf (i , • • • , «2n, ft) 

+ Pf (i , • • • , «2n, ft, c, d)Pf (io, ... , i 2n , a), (32) 
we obtain the relation between t^'* and Tn*: 

_ x s,t _ s+l,t-s,t+l _ „ _«,t+l-s+l,t I \ A*.* _s+l,t+l 

\ S+l,t+l-S,* _ + _ S+l,t-S,t+l , «,* *S + l,t+l /qq\ 

'V'n. 'n ~~ ^M'n 'n ^A' n 'n MiA 'n+l n— 1 ' V°°J 

where jz^a = z u -?a- Hence, from Theorem 2, we can easily arrive at the result. ■ 

Observing Theorem 3, we can regard the relation (29) as the discrete Lax pair whose wave 
function is the sequence of even-degree SOPs. Then we can directly obtain the following system. 

Theorem 4. From the compatibility condition of (29), we can derive the following nonlinear 
difference system: 

AS+l,t+l AS,t , r>s,t R s+l,t+l _ /-<S,i+l /-f«+l,t I r>S+l,t ns,t+l 

A n -°n+l D n ~ u n u n + U n U n , 

/(S + l,i/^S + l,t _ AS,tf*tS,t AS,t+l T)S,t + l _ AS,tT)S,t 

■^n u n-l — A n u n > ^n-l ~~ > 

W W = K+l^ = S&C^. (34) 

T/iis system has a special solution which can be written in terms of Pfaffians as (30). 

We can consider the system (34) as a discrete integrable system since this admits a discrete 
Lax pair by construction and we can put plenty arbitrary parameters into the solution (30). 

From the point of view of the r- function t^' 1 , the relation (24) which the elements of the 
Pfaffian satisfy is the same as the dispersion relation of the discrete coupled KP (dcKP) equa- 
tion [10] and indeed one can notice that the equation (33) is just the same as the dcKP equation. 
Under the choice of z as the spectral parameter, the dcKP equation yields the discrete Lax 
pair (29). Therefore the system (34) can be regarded as the nonlinear form of the dcKP. 
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4 An extension to 2 X 2 matrix form 

In the previous section, we propose the discrete integrable system in (2+l)-dimensional form (34) 
from the discrete spectral transformation of the SOPs. This system, by construction, is related 
only to the even-degree SOPs and not to odd-degree SOPs. Taking it into account that the 
integrable systems associated with both even and odd degree SOPs were already derived [1, 27], 
we can expect that there should exist a discrete integrable system in 2+1 dimension which is also 
connected to the discrete spectral transformations of both even and odd degree SOPs. In order 
to construct such system, as is described in the previous section, the ambiguity of odd-degree 
SOPs (3) is the barrier. Conversely, in this section, we shall construct such system not avoiding 
the ambiguity of odd-degree SOPs but utilizing this. 

Using the skew-inner product with the property (14), we introduce the new sequence 

of polynomials (z)}^ =0 : 



^2n( Z ) = "V, ^n+lW = ~^t, (35) 



where o , n ' t and a n ' are defined by 



*'* := Pf (0, . . . , 2n - 2, 2n) + (s/i + t\)r^, 



o* % := Pf (0, . . . , 2n - 1, In + l,z) + (s/x + t\)f%\ (36) 

respectively. Recall that r^'* and Tn* are defined by (22). here we remark that (p2n( z ) are 
no longer monic except for s = t = 0. Nevertheless, from Proposition 1, one can easily find 
{0n*(^)}n=o are also SOPs with respect to (-|-) s '*. 

As is seen in the previous section, in order to derive the discrete Lax pair of even-degree SOPs, 
that is the relation between Tn * and fn' , the equation (33) plays a central role. According to this 
observation, we wish to acquire the relation between Tn , fn , o~ S n and <5"«*. First we examine 
the properties of <7 n '* and a S n ' ■ 



Lemma 2. The function a n ' and a n ' defined in (36) hold the following properties: 
a s+i,t = pf ^ ... j 2n - 1, 2n + 1, /i) + (a^t + *A)r* +1 '* 



a s,t+l = pf ( 0) . . . , 2n - 1, 2n + 1, A) + (s/i + tX)r^ t+ \ 

a s+W = A -i P f ( 0) . . . f 2n, 2n + 2, /i, A) + ( Sf i + t\)r* +1 ' t+1 , 

<^ +1 '* = ^ lpf (°» • • • , 2n, 2n + 2, /z, z) + (sfi + t\)f s n +1 '\ 

°n +1 = z \ lpf (0, . . . , 2n, 2n + 2, A, z) + (s// + t\)f°> t+1 , 

K +W = A-^-^Pf (0, . . . , 2n + 1, 2n + 3, /x, A, z) + ( Sf i + a)f* +1 '* +1 . 

Proof. In a similar manner to the proof of Theorem 2, we can see 

Pf (0, . . . , k + 1, . . . , 2n + 1, fj,) = Pf (0, . . . , k, k + 2, k + 2, . . . , 2n, jj) 

= Pf (0, ...,k,k + 2, k + 2, ...,2n,n)- /i 2 Pf (0, ...,k,k, k + 2, ...,2n,n) 
= Pf (0, . . . fc, (Jfe + 2) - /u 2 • fc, k + 2, ...,2n,n) 
= Pf (0, 0, . . .k,(k + 2) - fi 2 ■ k, k + 2, ...,2n,(j,). 

Noticing (A; + 2) — fi 2 • k = k + 1 — fi ■ k, we get 

Pf(0, . . . , jfe + 1, . . . , 2n + 1, fi) = Pf (0, 0, ,2n,ii) - /xPf (0, 0, . . . , fc+1, ...,2n,fi) 
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= T s [Pf(0, . . . , k, . . . ,2n - 1)] - //T s [Pf(0, . . . , k + 1, . . . , 2n - 1)], 

where we have used the relation (26). Therefore with the symmetry fx A, we have 

T s [Pf (0, 1, . . . , 2n - 2, 2n)] = Pf (0, 1, . . . , 2n - 1, 2n + 1, /i) - 

T t [Pf (0, 1, . . . , 2n - 2, 2n)\ = Pf (0, 1, . . . , 2n - 1, 2n + 1, /x) - A<' m . 

In a similar manner, we can also obtain 

T s o r t [Pf (0, 1, . . . , 2n - 2, 2n)] = A^Pf (0, 1, . . . , 2n - 1, 2n + 1, /j) - (/x + A)t* +1 '* +1 , 
T s [Pf (0, 1, . . . , 2n - 1, 2n + 1, z)] = z^Pf (0, 1, . . . , 2n, 2n + 2, M , z) - 
T t [Pf(0, 1, . . . , 2n - 1, 2n + 1, z)] = x Pf(0, 1, . . . , 2ra, 2n + 2, A, 2) - Af*' m , 
r,or t rPf(0,l,...,2n-l,2n+l,Js)] 

= z-ipf (0, 1, . . . , 2ra + 1, 2n + 3, /i, A, z) - (/x + A)t£ +1 '' +1 . 

These results immediately lead us to the conclusion. ■ 

By virtue of Lemma 2, we can find the relation between r^' , Tn*, cr n ' and a" n ' • 

Theorem 5. The function r^'*, Tn *, On'* an d &rf defined in (22) and (36) satisfy the following 
four relations: 

\ s,t * s+l,t+l _ z S+l,t fiS,t+l _ S,t+l~S+l,t i \ .»+l,WiJ,t 

„ x _*,t *s+l,i+l _ _ s,t+l-s+l,t , \ s+l,*+l*s,t 

*ftA A fi' n +l w n-l ~ ^"n 'n 'n i A fJ.'n "n J 

^fj,,X^fj.o n+1 i n — -£\i n+ i o n ^M'n+l "n * A v°n 'n+1' 

„ \ _a,t *s+l,t+l _ 7l s+l,tis,t+l _„ s.t+l's+l,* , > s+l,*+l-s,t /o ? ^ 

Proof. Pfaffian identities (31), (32) and Lemma 2 directly lead us to the (37). This completes 
the proof. ■ 

If we take a = r, one finds that the equations (37) go to the discrete coupled KP equation (33). 
Therefore, these equations can be considered as the natural extension of the discrete coupled 
KP equation from the point of view of the SOPs. Hereafter we will call (33) as the extended 
discrete coupled KP (edcKP) equation. 

The edcKP equation gives us the contiguous relation among {(^^(z)}: 

Theorem 6. Let $ n ' (z) be the vector function given by 



LS,t I 

K 



The vector functions {&n (z)}^L satisfy the following two vector linear equations: 
z x ^\z) - z^^z) = z,, x A^ s n + _Y +1 (z) - BtfQtfiz), 

z„ tX $ s + l ' t+1 {z) ~ $n+i(z) = z x C s /^ n t+l (z) - z^ s /K +ltt , (38) 
where the coefficients A n '*, Bn' ! , C n ' , D n '* can be written explicitly in terms of Pfaffians: 
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B2' = (A-m) 



s,t s+l,t+l' 

s+l,t sl+l 
<?n 0~n 



s,t s+l,t+l 
Cn 0"ra 

y s+i,t s,t+i 



(a-m)- 1 



+l,Li,t+l 



T n+1 T ' 



i.t s+l,t+l 



= (A- M r 1 



5+l,i s,t+l 

. s,t s+l,t+l 
\ T n+l T n 



s,t+l s+l,t 

■ s,t s+l,t+l 
\ T n+l Tn 



s,t+l s+l,t ' 
T n+1 T n 
s,t s+l,t+l 



(39) 



/ 



Proof. The definitions of }™ =0 (35) and the edcKP equation (37) directly lead us to the 
conclusion. This completes the proof. ■ 

We can regard (38) as the discrete vector Lax pair of SOPs of both even and odd degree. 
Therefore the discrete Lax pair shows us the following systems associated with SOPs of both 
even and odd degree, which can be regarded as a discrete integrable system along with (34): 

Theorem 7. From the compatibility condition of the discrete vector Lax pair (38), we can derive 
the nonlinear system in 2 + 1 dimension which takes the 2x2 matrix form: 



+l,tpis+l,t _ \8,t(iB+l,t a s,t+l-pjS,t+l _ is.ins,! 



P>S+l,t-pjS+l,t _ p>S,t pvS,* 

D n — ■ D n +1- L 'n ) 



Bs,*+1-T»s,t+1 _ r,s,t n si 



(40) 



Special solutions of the nonlinear system (40) can fre written in terms of Pfaffians given by (39) . 

One can easily find that this system is the 2x2 matrix extension of the system (40) and 
this is the nonlinear form of the edcKP equation. This system is clearly connected to the SOPs 
of both even and odd degree and from its construction, we can regard the system (40) as the 
discrete coutnerpart of the Pfaff lattice in (2 + l)-dimensional form. 



5 Concluding remarks 

In this paper, we proposed the discrete spectral transformation of SOPs. From the discrete 
spectral transformation of SOPs, we derived the discrete integrable systems both in (1 + 1)- 
dimensional form and in (2 + l)-dimensional form. In (1 + l)-dimensional case, the system gives 
the integrable discretization of the Pfaff lattice, though the system is of course not local. In 
(2 + l)-dimensional case, the system is local and corresponds only to the even-degree SOPs. This 
system is also shown to be connected to discrete coupled KP equation. The discrete coupled 
KP equation is known to be equivalent to the discretized Hirota-Ohta equation [11] in the 
bilinear form levels. Recently, the discrete integrable system is proposed from the auxiliary 
linear problems for the Hirota-Ohta (or DKP) system [4, 24]. 

We further extended the discrete integrable system associated with even-degree SOPs to 
the system associated with both even and odd degree SOPs. This system can be described in 
2x2 matrix form, which is related to an extension of the discrete coupled KP equation. This 
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extension has never been considered before, although it is quite natural from the aspects of 
SOPs. 

Finally, let us mention for the relationship between our results and the random matrix theory. 
As a by-product of the discrete transformation of SOPs, we found the compact factorization of 
the kernel function of the correlation function for the orthogonal or symplectic ensembles in 
terms of the two even-degree SOPs with respect to the different weights. This factorization is 
not known as far as the authors know and indeed might be a new one. Moreover, it is known 
that a matrix integral over orthogonal and symplectic ensembles of Hermitian matrices can be 
considered as a r-function of the coupled KP hierarchy [12], while our discrete integrable sys- 
tems (34) and (40) are directly connected to discrete coupled KP and its extension, respectively. 
We hope these results will bring a deeper understanding of the random matrix theory and several 
applications in its works in the near future. 
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